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The oscillator strengths of the H^" molecular ion, lsa g — 2pa u , lsa g — 2pn u are calculated within the 
Born-Oppenheimer approximation. The variational expansion with randomly chosen exponents has 
been used for numerical studies. The oscillator strengths obtained for the transitions lsa g — 2pa u , 
lsag — 2pir u of Hj are accurate up to ten significant digits. Results are given for the internuclear 
distances between 0.10 and 20.0 a.u. 



I. INTRODUCTION 



Q_i, Many theoretical calculations on the oscillator strengths of the Hj molecular ion have been performed since the 

J_i ■ publication of the first [1] of Mulliken's series of papers on the subject. Bates and his co-workers had done a great deal 
of work on this molecular system, with special emphasis on the electronic transition probabilities for fixed internuclear 
separation R for a large number of electronic transitions [2|-|6[ • The problem was formulated in terms of Jaffe's solution 
[7[ and led to certain integrals in the expression for the transition probabilities which required numerical evaluation. 

^ . Moreover, using the Hylleraas solution [8| the matrix elements of the electric dipole moment for transition between 
bound electronic states of the Hj" molecular ion for any fixed internuclear separation R have been obtained in analytical 
form 9]. In the last years, with the aid of the new types of variational expansions based on the randomly chosen 

^V exponents for the Coulomb three-body problem [lCt lll| , within the Born-Oppenheimer approximation, very accurate 
f^ non-relativistic energies and the improved relativistic corrections of ma 4 and ma 6 orders for the ground state of the 

■^H ' ion, as well the more accurate static dipole polarizability for the lsa g electronic state of Hj molecular ion have been 
obtained in [J^Wli}. 

In this work our goal is to obtain more accurate values for the oscillator strengths of the lsag — 2pa u , ls<j g — 2p-K u 



transitions of HJ" molecular ion within the Born-Oppernheimer approximation than those in [2j and 15 1 



We want to show that the use of the variational expansion suggested in this work allows an analytical evaluation 
of the matrix elements and can provide us with very accurate data. In future studies it can be used to obtain more 
^_i_' . accurate values of moments functions Sk{R) as a function of internuclear separation R for the lsa g electronic state 
of the HJ" molecular ion in the Born-Oppenheimer approximation. 

o 

, ■ II. THEORY 



In what follows we consider the Coulomb two-center problem with the nonrelativistic Hamiltonian 



- '-I 

C^ 2m r\ T2 

where r\ and r-i are the distances from an electron to nuclei 1 and 2, respectively. 
The oscillator strength fo n (details can be found in [l|, 0]) can be expressed as 

hn = \G\Q\ 2 {E n -E Q ) (2) 

where G is the orbital degeneracy factor, and G — 1 for a states, G = 2 for it states. Here Q is the transition moments 
and can be written in the form 

Q = (0\r\n) , (3) 
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where r is a distance of the electron from the center along the internuclear axis: 

1 



., , V 2r l r 2 + 2r l R2 + 2r 2 R2 ~ r l - r 2 - RA ■ 

III. VARIATIONAL APPROXIMATION AND NUMERICAL RESULTS 

In order to get a precise solution for the Schrodinger equation 

»o(r) = Eo*o(r) (4) 






2m 

we use the variational approach. A variational wave function for the electronic states of the H 2 should be symmetrized 
and is constructed as follows: 

oo 

tf (ri, r 2 ) = e ™*Y |m| J^ C 4 (e~ ai, ' 1 ~ ftr2 ± e -^ ri -°" r2 ), (5) 

where (+) is used for the gerade electronic state and (— ) is for the ungerade electronic state, respectively. Parameters 
a, and /% are generated in quasi-random manner [10L lllfl 

a l =[-i{i + l) y /p^\(A 2 -A 1 ) + A 1 (6) 

\_x\ designates the fractional part of x, p a is a prime number, an interval [^4i, A2] is a real variational interval, which 
has to be optimized. Parameters Pi are obtained in a similar way. 

Numerical evaluation of the matrix elements for operators in ([2]) and ([3]) is expounded in the Appendix. 

In Tables I the energy differences and and oscillator strengths as a function of internuclear separation R for 
lser ff — 2pa u transitions (E\ — E and / i) and lsa g — 2pn u transitions (_E 2 — E and /02) and comparisons with 
the previous results are presented. The estimated accuracy of obtained values for the transitions is ten significant 
digits. In Table II the convergence of the lsa g — 2pa u and lsa g — 2pit u transitions with an increase in the expansion 
length of the \E r (equation ([5])) is demonstrated for some values of internuclear distance R. 

In order to get accurate results we use three sets of basic functions of the type ([5|) (in a spirit of Ref.[l(|) for small 
values of internuclear distance R, two sets for intermediate and large values of R, respectively. Total number of the 
basic function varies from N = 90 to N = 180. In Table II in the second and third columns the numbers of the basis 
functions employed (equation (|5|)) for the electronic states 2pa u and 2p , n u of HJ" are presented, respectively. In our 
calculations arithmetics of sextuple precision (about 48 decimal digits) implemented as a FORTRAN90 module has 
been used. 

In all tables the factor x in the brackets means 1Q X . Atomic units are used throughout. 

IV. CONCLUSION 

The oscillator strengths of the H^ molecular ion, H^, lsa g — 2pa u , lsa g — 2pir u have been calculated accurately 
within the Born-Oppenheimer approximation. The variational expansion with randomly chosen exponents has been 
used for numerical studies. This type of expansion allows us to use few number of basic functions, and even the non- 
relativistic energy for the lsa g electronic state of the ion can be obtained to 10 _20 a.u. accuracy when the number of 
basis functions is less than N — 100 |14l |. 

Previous calculations for the lscr g — 2pa u transitions of the ion for a fairly-restricted internuclear distance R are 
by Mulliken [lj, Bates [2], and Dalgarno and Poots |h|. Mulliken [l] used the LCAO approximation and Bates [2] 
employed the exact two-center wave function for the numerical integral evaluation. Dalgarno and Poots |6j employed 
the variational method in the calculation. The calculations for the lsa g — 2pir M transition of Hjf are by Bates, 



Darling, Hawe and Stewart (4J, Herman and Wallis [9J and Ramaker and Peek [TaL Bates et al [J] used also the 
exact two-center wave function for the numerical integration. Herman and Wallis [9j obtained the matrix elements 
for the transitions between bound electronic states of the ion for any fixed internuclear distance R in an analytical 
form using the Hylleraas solution for the wave functions. Ramaker and Peek [15| performed the calculations for the 
dipole strengths of the Hj molecular ion as functions of wide range of internuclear separation R for all transitions 
the states which correlate to a proton-plus- hydrogen- atom and have principal quantum number n = 1,2 or 3, using 
the exact electronic two-center wave function. 



In Table III, we compare our results with earlier ones, which demonstrate the superiority of the newly obtained 

■ g - 2pa ui lscr g - 2pir u of H+ 



results. The results obtained for the transitions ls<7„ — 2pa u , lser„ — 2pir u of HJ are accurate up to ten significant 



digits, that is, six digits more than in 0, LL5[. With the help of the results obtained in this Letter, one can recalculate 
the mean excitation energy for the lsa g electronic state of the Hj" molecular ion employing the functional form of the 
moment functions Sk(R) [16l4l8| |. This work is in progress now. 

Appendix: Analytical evaluation of the matrix elements 

The calculation of the matrix elements is reduced to evaluation of integrals of the type 

T lm (a,/3)= I ' r[- l r™- l e- a ^-^d?v. (A-l) 



Integers (I, m) are, in general, non- negative, but in case of singular matrix elements one of the indices can be negative. 
The function r o can be easily obtained 

4-ir p~P R — p- aR 
r 00 (a 1 ^R) = - a2 _ p2 , (A-2) 

where R is the distance between nuclei, then Fi m (a, /3; R) for non-negative (l,m) may be generated from (IA-21) by 
means of relation 

d V / d 



r im (a,/3;i?)= (-— 1 (-— ) r o(a,/3,£). (A-3) 



Integral r_ 10 (a, /?; R) is expressed by 



27T ( 

r_ 1)0 (a,i8; R) = — |e /3i? E 1 (i 1 '(a + /3)) + e^ln^a + /3) 
-eP R [E x (R(a - /?)) + lnJJ(a - 0)] }. 



(A-4) 



Worthy to note that a function in square brackets is analytic when argument is zero. Integrals T_i !m are generated 
from r_io similar to (|A-3|) : 

T_ hm (a,l3;R)= (- — J r_i, (o,j8,U). (A-5) 

Function Ei(z) encountered in (|A-4=|) is the exponential integral function [19j : 

Ex(z) =r(0,z)= / t^e^dt. 
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TABLE I: Data on the lsa g — 2pa u , lsa g — 2pir u transitions of H^ 







Bates 


2] 




Ramaker and Peek [15] 




R 


Ei — Eq 


/oi 


/oi 


E-2 — Eo 


/02 


/02 


0.10 


0.147 757 4626(1) 




0.1419461197 


0.147 857 3935(1) 




0.281958 4589 


0.20 


0.142 594 2876(1) 


0.150 


0.149 9661487 


0.142 9931822(1) 




0.292 176 8241 


0.30 


0.136 065 8191(1) 




0.161227 2533 


0.136 9601755(1) 




0.305 026 4125 


0.40 


0.128 996 9831(1) 


0.175 


0.174 800 4051 


0.130 578 7951(1) 




0.318 8716774 


0.50 


0.121810 2535(1) 




0.190 006 5537 


0.124 264 4532(1) 




0.332 827 2240 


0.60 


0.114 717 4343(1) 


0.207 


0.206 235 8257 


0.118 218 5058(1) 




0.346 408 7672 


0.70 


0.107 822 0761(1) 




0.222 877 7173 


0.112 529 8372(1) 




0.359 356 2504 


0.80 


0.101173 4174(1) 


0.240 


0.239 317 2609 


0.107 228 2899(1) 




0.371 539 1200 


0.90 


0.947 946 4449 




0.254 970 6910 


0.102 313 0356(1) 




0.382 902 9462 


1.00 


0.886 972 6883 


0.269 


0.269 336 9702 


0.977 678 3571 


0.393 


0.393 4381322 


1.30 


0.721475 3157 




0.301 731 0394 


0.861082 3079 




0.420 300 6176 


1.50 


0.625 8216700 


0.314 


0.313 859 3007 


0.797 263 3391 


0.435 


0.434 630 6848 


1.80 


0.503 5781216 


0.321 


0.320 620 5783 


0.717 917 7718 


0.451 


0.4515518442 


2.00 


0.435 099 8223 


0.319 


0.319 763 3919 


0.673 862 3946 


0.460 


0.460 186 9855 


2.30 


0.349 292 2274 




0.313 4811372 


0.618 286 6408 




0.469 722 4456 


2.50 


0.3017514622 


0.307 


0.3071116776 


0.587 007 7894 


0.474 


0.474 048 0566 


2.80 


0.242 384 3555 


0.297 


0.295 4491732 


0.5471114888 


0.478 


0.477 8115429 


3.00 


0.209 477 8640 


0.289 


0.286 648 3110 


0.524 452 3451 


0.479 


0.478 648 2112 


3.50 


0.145 359 8123 


0.264 


0.262 040 2032 


0.478 772 9801 


0.475 


0.475 449 5457 


4.00 


0.100 534 2444 


0.238 


0.234 600 2369 


0.445 260 2396 


0.465 


0.465 523 7278 


4.50 


0.691 110 5645(-l) 


0.295 


0.205 243 3086 


0.420 759 4992 


0.450 


0.449 9718773 


5.00 


0.471 286 8194(-1) 


0.175 


0.175 094 9934 


0.403 035 4828 


0.430 


0.430144 3360 


5.50 


0.318 483 6379(-l) 


0.144 


0.145 512 8632 


0.390 423 8948 


0.408 


0.407 637 0403 


6.00 


0.213 251 5610(-1) 


0.116 


0.117 844 6702 


0.381 644 1528 


0.384 


0.384106 9351 


6.50 


0.141 551 5325(-l) 


0.090 


0.931 480 9973(-l) 


0.375 699 5692 


0.361 


0.361 026 7789 


7.00 


0.932 229 1695(-2) 


0.069 


0.720 267 2925(-l) 


0.371819 5628 


0.339 


0.339 507 0520 


7.50 


0.609 732 1052(-2) 


0.090 


0.546 252 3165(-1) 


0.369 418 4230 


0.320 


0.320 242 4604 


8.00 


0.396 437 2998(-2) 


0.039 


0.407 352 9645(-l) 


0.368 059 7573 


0.304 


0.303 561 4626 


8.50 


0.256 445 4385(-2) 


0.028 


0.299 374 4120(-1) 


0.367 424 2386 


0.289 


0.289 523 5398 


9.00 


0.165 162 3429(-2) 


0.020 


0.217256 4226(-l) 


0.367 2814212 


0.278 


0.278 018 0087 


9.50 


0.105 969 2322(-2) 


0.016 


0.155 938 7742(-l) 


0.367 466 5422 


0.269 


0.268 8416239 


10.00 


0.677 660 3412(-3) 


0.011 


0.110 852 7572(-l) 


0.367 862 4388 


0.262 


0.261750 5281 


11.00 


0.274 807 5098(-3) 




0.546 544 7478(-2) 


0.368 978 5760 




0.252 8151953 


12.00 


0.110 412 3031(-3) 




0.262 186 7652(-2) 


0.370 217 2763 




0.249 293 6287 


13.00 


0.440 274 9253(-4) 




0.122 982 9099(-2) 


0.371377 5187 




0.249 503 7766 


14.00 


0.174 459 4699(-4) 




0.566 128 7073(-3) 


0.372 370 6777 




0.252 0276563 


15.00 


0.687 612 7753(-5) 




0.256 469 7079(-3) 


0.373 170 2448 




0.255 725 8261 


16.00 


0.269 771 2233(-5) 




0.114 593 6419(-3) 


0.373 783 9614 




0.259 757 3081 


17.00 


0.105 416 2788(-5) 




0.505 885 9722(-4) 


0.374 236 6462 




0.263 577 2171 


18.00 


0.410 477 9229(-6) 




0.220 971 2567(-4) 


0.374 5591664 




0.266 896 4478 


19.00 


0.159 336 3912(-6) 




0.956 149 8035(-5) 


0.374 781 7742 




0.269 614 3011 


20.00 


0.616 776 4338(-7) 




0.410 256 5319(-5) 


0.374 930 7098 




0.271746 9205 



TABLE II: Convergence of the \so g — 2pa u and lsa g — 2pn u transitions with an increase in the basic-set size of \& 



R(a-u-) N foi fm 

0.10 110 65 0.141946119 692 60 0.281958 458 927 66 

125 80 0.141946119 692 60 0.281958 458 927 67 

140 95 0.141946119 692 60 0.281958 458 927 67 

2.00 90 68 0.319 763 391895 63 0.460186 985 489 56 

120 78 0.319 763 391895 63 0.460186 985 489 56 

140 98 0.319 763 391895 63 0.460186 985 489 56 

20.0 90 105 0.410 256 531 872 29(-5) 0.271746 920 513 63 

110 125 0.410 256 531 861 97(-5) 0.271746 920 515 47 

130 160 0.410 256 531 861 96(-5) 0.271746 920 515 47 



TABLE III: Comparison with earlier calculations at a bond length R = 2.0 a.u. 





/oi 


/02 


Mulliken [1] 


0.31 




Dalgarno and Poots [6] 


0.319 




Bates et al [2, 4] 


0.319 


0.460 


Ramaker and Peek [15] 


0.319 


0.460 


This work 


0.319 763 3919 


0.460186 9855 



